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Abstract 

The problem of a test body in the Schwarzschild geometry is investigated in a Keplerian limit. 
Beginning with the Schwarzschild metric, a solution to the limited case of approximately elliptical 
(Keplerian) motion is derived in terms of trigonometric functions. This solution is similar in form 
to that derived from Newtonian mechanics, and includes first-order corrections describing three 
effects due to general relativity: precession; reduced radial coordinate; and increased eccentricity. 
The quantitative prediction of increased eccentricity may provide an additional observational test 
of general relativity. By analogy with Keplerian orbits, approximate orbital energy parameters 
are defined in terms of a relativistic eccentricity, providing first-order corrections to Newtonian 
energies for elliptical orbits. The first-order relativistic equation of orbit is demonstrated to be a 
limiting case of a very accurate self-consistent solution. This self-consistent solution is supported 
by exact numerical solutions to the Schwarzschild geometry, displaying remarkable agreement. A 
more detailed energy parameterization is investigated using the relativistic eccentricity together 
with the apsides derived from the relativistic effective potential in support of the approximate 
energy parameters defined using only first-order corrections. The methods and approximations 
describing this Keplerian limit are applied to more general static spherically-symmetric geometries. 
Specifically, equations of orbit and energy parameters are also derived in this Keplerian limit for 
the Reissner-Nordstrdm and Schwarzschild-de Sitter metrics. 
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I. INTRODUCTION 



The problem of a test body in the Schwarzschild geometry is of fundamental importance 
in understanding orbital characteristics due to general relativity. Detailed analyses and ap- 
proximate analytical solutions of this geometry exist [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 
14, 15, 16, 17, 18], emphasizing a first order approximate rate of precession of perihelia of 
Mercury and the unstable circular orbit and relativistic capture representative of extreme 
astrophysical environments. The general problem has been solved in terms of WeierstraB 
elliptic functions by Hagihara [19], and more recently by Kraniotis and Whitehouse [20]. 
Ashby [21] has solved the limited case of approximately Keplerian motion in terms of Jaco- 
bian elliptic functions, "...which are closer in spirit to trigonometric functions, with which 
most are familiar." In this paper, beginning with the Schwarzschild geometry, a solution 
to the limited case of approximately Keplerian motion is derived in terms of trigonomet- 
ric functions. This approximate solution lends itself to easy comparison with the familiar 
Keplerian orbits (ellipses) of Newtonian mechanics and is consistent with both the reduced 
radius of circular orbit, commonly derived from the Schwarzschild effective potential, and 
the observed rate of precession of perihelia of the inner planets, commonly derived pertur- 
batively. Additional insights regarding the sizes and shapes of bound relativistic orbits are 
provided, including a relativistic correction to eccentricity which may be subjected to obser- 
vational tests. By analogy with Keplerian orbits, a relativistic eccentricity is used to define a 
Schwarzschild energy parameter, providing corrections to Newtonian energies for Keplerian 
orbits. The relativistic equation of orbit together with the energy parameterization com- 
prise a simple model that is useful for a qualitative and quantitative understanding of the 
corrections to Keplerian orbits due to general relativity. This model is easily extended to 
include more general static spherically-symmetric geometries. Specifically, models are also 
derived for the Reissner- Nordstrom and Schwarzschild-de Sitter metrics. 

The methods and approximations describing a Keplerian limit are detailed in Sec. II. 
Beginning with the Schwarzschild geometry, an approximate equation of orbit is derived 
in terms of trigonometric functions. When compared to that of a corresponding Keplerian 
orbit, this equation of orbit clearly displays three characteristics of relativistic orbits: pre- 
cession; reduced radial coordinate; and increased eccentricity. These characteristics arise 
as first-order relativistic corrections to the familiar equation of orbit describing Keplerian 
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orbits. This solution is found to be valid for near-circular orbits requiring only small rela- 
tivistic corrections. Predictions of relativistic precession and reduced radius of circular orbit 
are in agreement with known results. This provides confidence in a new relativistic correc- 
tion to eccentricity. (The possibility of observing relativistic corrections to eccentricity is 
discussed briefly in Sec. II A.) In addition, first-order relativistic corrections to Keplerian 
apsides are predicted, resulting in the conclusion that the overall size of a relativistic orbit 
is smaller than a corresponding Keplerian orbit. By analogy with Newtonian mechanics, 
Schwarzschild energy parameters are defined using the virial theorem for circular orbits, 
and using a relativistic eccentricity for noncircular orbits, resulting in first-order relativistic 
corrections to Newtonian energies for Keplerian orbits. This simple model is substantiated 
by a more detailed investigation in Sec. III. 

Again beginning with the Schwarzschild geometry, the derivation of a very accurate self- 
consistent relativistic equation of orbit is detailed in Sec. III. This self-consistent equation 
of orbit is identical in form to the more approximate equation of orbit derived in Sec. II 
and predicts the same orbital characteristics. However, corrections to Keplerian orbits due 
to general relativity are described more accurately. For example, the predicted radius of 
circular orbit is identical to that derived by minimizing the relativistic effective potential. 
This solution is accurate in predicting long-term behavior of Schwarzschild orbits for a 
large parameter space, as is demonstrated by comparisons with exact numerical solutions. 
The more approximate equation of orbit derived in Sec. II is shown to be a limiting case 
of this self-consistent solution. A more detailed energy parameterization is investigated for 
comparison with the simple parameterization of Sec. II. This parameterization is constructed 
using the relativistic eccentricity derived from the self-consistent equation of orbit together 
with the relativistic apsides derived from the Schwarzschild effective potential. The energy 
parameters of Sec. II approximate well the results from this more detailed parameterization, 
lending value to the simpler approach and more approximate results. 

The methods and approximations describing a Keplerian limit to the Schwarzschild ge- 
ometry are applied to a more general class of static spherically-symmetric geometries in 
curvature coordinates [4, 5, 6, 22, 23, 24]. Specifically, the path of a small test mass in a 
static spherical spacetime is taken to be described by the metric 
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where di?^ = dO'^ + sin^^df/)^, and 

r r"' 3 

The Reissner-Nordstrom (A = 0) geometry is considered in Sec. IV, followed by the 
Schwarzschild-de Sitter (Q = 0, A > 0) geometry in Sec. V. These geometries have been 
shown to meet all of the conditions for physical acceptability [25, 26, 27, 28, 29]. In each 
of these cases the relativistic correction due to matter is considered to be the dominant 
contribution. The resulting orbital equations are of the same form as that derived for the 
Schwarzschild geometry, including additional corrections for each geometry. For examples: 
there is an additional contribution to relativistic precession due to charge opposite in direc- 
tion to the contribution due to matter; there is also an additional contribution to relativistic 
precession due to A, but in the same direction as the contribution due to matter. Nu- 
merical studies are provided only for the Schwarzschild geometry in order to establish the 
self-consistent approach. A self-consistent treatment of the Schwarzschild-de Sitter geometry 
is intractable and is not pursued. A concise summary of the results for the three geometries 
is given in Sec. VI. 



II. SCHWARZSCHILD ORBITS IN KEPLERIAN LIMIT 

The path of a small test mass near a spherically-symmetric central mass M is uniquely 
described by the Schwarzschild geometry [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 30, 31]. 

ds' = (1 - r.^McMt' - - rM^2^ (3) 

1 - r^j/r 

where di7^ = d9^ + sin^Odip"^. The singularity at the Schwarzschild radius r^^ = 2GM/c^ 
is irrelevant in the present context in which a solution far from the central mass is sought. 
Consider orbits in the plane defined by ^ = 7r/2. Parameterize timelike geodesies with 
ds^ = c^dr^, where r is the proper time along the path of a test particle. Then, with 
= dx'^/dr, the equations of motion may be expressed as 

e = r^vb, (4) 
fc^(l-Wr)t, (5) 
i^e-l)c' = lr' + V,s, (6) 



4 



where a Schwarzschild effective potential is defined as 

rr^K,--^+i;i-%^. (7) 

The parameter = i'^/GM is the radius of a circular orbit for a nonrelativistic particle 
with the same angular momentum i, and the mass-related relativistic correction parameter 
is defined as 

_ [UIVl ' ^ 



/GMV ^ 1 

; 2 ■ 



(8) 



The Newtonian effective potential is recovered in the limit e^^ — * 0, 



r'j-'V^s^-'- + ^^. (9) 

Time is eliminated from Eq. (6) using r = —Mu/d(p, where u = 1/r. An equation for the 
trajectory of a test particle is then obtained by differentiating once more with respect to ip, 

/^- + - = 1 + 36m(-)'. (10) 

dip r r \ r / 

The conic-sections of Newtonian mechanics are recovered by setting = in Eq. (10), 

-r^ \ = 1 =^ — = l + ecosv2, (11) 

d(/? r r r 

where e is the eccentricity of the orbit and is taken to be positive or zero. 

For approximately Keplerian orbits it is convenient to linearize the equation of motion 

(10) by making the change of variable 

^^r 1 / X 

- - 1 = - < 1, 12 
r a 

so that the last term on the right-hand-side of Eq. (10) may be approximated as 

(13) 



rc\2 . 2 



a 



Equation (10) may now be expressed as 



1 d^ 1 1-66^,1 



1. (14) 



3% d(p'^ a Sejv/ a 

Defining = (1 — 6e^^)/3e^j and making the additional change of variable a = (1 — 6ej^j)^(p 
results in the familiar form 

+ ^ - 1. (15) 
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The solution is described by Eq. (11) 



a 



1 + A cos a, (16) 



where A is an arbitrary constant of integration. In terms of the original coordinates, Eq. (16) 
becomes 

T 

— ~ 1 + e COS /to9, (17) 
r 



where 



r.\=^, (18) 



1 - 3% ■ 



(19) 



K = (l-6e^,)5. (20) 

According to the correspondence principle, the solution must reduce to that of Newtonian 
mechanics in the limit ej^.^ 0. Therefore, 

3eMv4 = e (21) 

is identified as the eccentricity of Newtonian mechanics (11). Then, to first order in e^, 
Schwarzschild orbits in this Keplerian limit are described by Eq. (17), where 

fe^r,(l-3eA,), (22) 
e^e(l + 3ejv,), (23) 
~K^l-3ej,j. (24) 

Therefore, Schwarzschild orbits in this Keplerian limit may be expressed concisely as 

^l + e{l + 3e^j) cos (1 - Se^)^. (25) 

A systematic verification may be carried out by substituting (25) into (10), keeping terms 
of orders e, e^^.^, and eej^j only. However, the justification for discarding the term nonlinear in 
eccentricity is the correspondence principle. Arguments concerning which terms to discard 
based only on direct comparisons of relative magnitudes of higher-order and lower-order 
terms lead to contradictions. Rather, the domain of validity is expressed by subjecting the 
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solution (25) to condition (12) for the smallest value of r. Evaluating the equation of orbit 
(25) at pericenter (r = f_) results in 

rc _ 1 + e(l + Se^) . . 

so that, according to (12), the domain of validity is given by 

e(l + 36^,) + 2(3e^,) « 1. (27) 

Therefore, the relativistic eccentricity e = e(l + Sej^j) <^ 1; the Schwarzschild equation of 
orbit (25) is limited to describing relativistic corrections to near-circular, elliptical orbits. 
Also, the relativistic correction 6e^j ^ 1; the Schwarzschild equation of orbit (25) is valid 
only for small relativistic corrections. 

Although the more general solution, Eq. (17) together with (18)-(20), contains higher- 
order terms in e^^j, it is only consistent to first order. A self-consistent solution that is 
consistent to all orders in e^^^ and accurately predicts long-term orbital behavior is presented 
in Sec. III. However, the first-order relativistic equation of orbit (25) is similar in form 
to that describing Keplerian orbits (11) and is useful for a qualitative and quantitative 
understanding of relativistic corrections to Keplerian orbits. When compared to Keplerian 
orbits, Schwarzschild orbits clearly display three characteristics: precession; reduced radial 
coordinate; and increased eccentricity. 



A. Characteristics of Schwarzschild Orbits 



The approximate equation of orbit (25) predicts a shift in apside through an angle 

= 2n{k^^ - 1) (28) 
^ 3eM(2vr) (29) 

per revolution. This first-order prediction is in agreement with existing perturbative calcu- 
lations [4, 22, 32], as well as with other calculations [14, 15, 16, 17, 18], and is well-known to 
be in agreement with the observed precession of perihelia of the inner planets [8, 13, 33, 34]. 
Precession due to relativity is illustrated in Fig. 1. The rate of precession is exaggerated 
by the choice of relativistic correction parameter (e^^ = 0.1) for purposes of illustration. 
However, relativistic orbits precess for smaller (non-zero), reasonably chosen values of as 
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well; Figure 1 correctly illustrates that relativistic orbits precess, a characteristic which is 
not present in Keplerian (central-mass) orbits. 

The approximate equation of orbit (25) predicts a reduced radius of circular orbit, 

5fc = fc-rc^ -Se^Tc. (30) 

This prediction is in agreement [1, 7, 35] to first order with that obtained by minimizing the 
Schwarzschild effective potential (7), 

Rc = lr, + lr,Vl-l2e,,. (31) 

(A distinction is made here between the radius of circular orbit as determined from the 
Schwarzschild effective potential -Rc and that determined from the equation of orbit fc.) For 
ejyj <^ 1/12, the radius of circular orbit is predicted to be reduced, 

6Rc = Rc-rc^ -3ejv^rc. (32) 

The Schwarzschild effective potential (7) is compared to that derived from Newtonian me- 
chanics (9) in Fig. 2. 

The approximate equation of orbit (25) provides the further insight that a non-circular 
Schwarzschild orbit is predicted to be smaller than a corresponding Keplerian orbit for a 
large range of eccentricities. The relativistic apsides may be expressed approximately as 

f±/r± ^ 1 - 3e,,-^, (33) 

where r± = rc/(l =F e) denotes the Keplerian apocenter (-I-) and pericenter (— ). The rela- 
tivistic apocenter is equal to the corresponding Keplerian apocenter for e+ ~ 1/2. For all 
nonzero values of e^^j, the apocenter is reduced for e < e+ and increased for e > e+. The 
relativistic pericenter is reduced for all nonzero values of e^/. Therefore, for e < e+, the 
radial coordinate is reduced; the observed overall size of a relativistic orbit is smaller than 
a corresponding Keplerian orbit for e < e+. Using a perturbative treatment, Nobili and 
Roxburgh identified a correction of the same order of magnitude [14] : 

"There is a constant part, causing a variation in the average size of the classical 
orbit of the order of 

V « /c c^a 
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e.g. ~ 6 X 10 ^ for Jupiter and smaller for the outer planets; it is anyway smaller 
than the present accuracy in a and therefore can be neglected." 

(Therein, a refers to the Keplerian semimajor axis.) It is possible that future experiments 
will be sensitive to this first-order relativistic correction to the size of orbits. 

These characteristics are illustrated in Fig. 3, in which orbits derived from the 
Schwarzschild geometry (25) and Newtonian mechanics (11) are compared. These char- 
acteristics are exaggerated by the choice of parameters (e = 0.4, ej,^ = 0.1) for purposes of 
illustration, and precession has been removed {k 1) in order to emphasize the size and 
shape of the Schwarzschild orbit. However, Schwarzschild orbits have reduced radial coordi- 
nates for smaller (non-zero), reasonably chosen values of e and ej^j as well; Figure 3 correctly 
illustrates that Schwarzschild orbits are smaller than corresponding Keplerian orbits. 

When compared to Keplerian orbits, the Schwarzschild geometry predicts orbits that are 
more eccentric, 

5e = e — e ~ 3e^e. (34) 

This characteristic is also displayed in Fig. 3, in which it is noticeable that the semi-minor 
axis is reduced more than the semi-major axis. This quantitative prediction may be sub- 
jected to observational tests. It is possible that part of the discrepency between observed 
and calculated eccentricities in many astrophysical systems may be accounted for by this 
contribution. Bosch et al. [36] find galaxy and cluster substructure eccentricity distributions 
to be strongly skewed toward high eccentricity. Also, the discrepancy in the eccentricity of 
the massive binary blackhole system in OJ 287 is found [37] to be of the order expressed 
by Eq. (34). Champion et al. [38] have identified an eccentric binary millisecond pulsar 
(PSR J1903+0327) "...that requires a different formation mechanism..." in order to recon- 
cile the large eccentricity and short spin period. 

B. Schwarzschild Energy Parameters 

For a particular Keplerian orbit, identified by total energy, it is necessary to identify the 
corresponding Schwarzschild energy for an orbit described by the same angular momentum. 
This provides a relativistic correction to Newtonian energies for Keplerian orbits. This is 
also useful when comparing orbital properties using energy diagrams. The generalized virial 
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theorem [39, 40, 41, 42] provides this relation for circular Schwarzschild orbits, and is easily 
extended to more general metrics. Referring to (6) and (7), a Schwarzschild potential energy 
parameter for a circular orbit is defined as 

~ GM GMrl 

' c ' c 

where fc is the radius of a relativistic circular orbit (22), and e^^ is defined in Eq. (8). 
According to the virial theorem, 

~ GM GMrl 



where Tc is a Schwarzschild kinetic energy parameter. Therefore, using Ec = Tc + 



''^^'^-.J) (38) 



^c(l + 2e^,), (39) 



where Eq. (22) is used in the last step, and = —GM/2rc is the total energy per unit mass 
for a circular Keplerian orbit of radius Tc- This may instead be expressed as 

E,~2ej,j\E,\. (40) 

Therefore, 6Ec ~ — 2e^.j|£'c| < 0; The energy of a circular Schwarzschild orbit (e = e = 0) is 
less than that of a corresponding Keplerian orbit. This relation is displayed in Fig. 4. 

More generally, the quantitative prediction of increased eccentricity provides an avenue 
by which Keplerian and Schwarzschild energies may be compared. The total energy per unit 
mass for a Keplerian orbit may be expressed in terms of the eccentricity as [43, 44] 

E/E, = l~e\ (41) 

Analogously, an approximate Schwarzschild energy parameter is defined by ansatz to be 

E/E^ ^ 1 - (42) 

The Schwarzschild energy parameter (42) is then expressed in terms of the Keplerian energy 
using Eqs. (23), (39) and (41), 

E ^ E - 2e,i{l - A^)\E,\. (43) 
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(This ansatz is investigated in Sec. IIIB, wherein (43) is shown to approximate well a more 
detailed parameterization and is a good approximation for near-circular orbits.) Therefore, 
6E ^ ~2ej^j\Ec\{l — 4:6^) < 0; The Schwarzschild energy is less than the corresponding 
Keplerian energy until the energies become equal for cq ~ |. The Schwarzschild energy 
parameter becomes greater than the corresponding Keplerian energy for e > cq. For near- 
circular orbits the Schwarzschild energy parameter always lies below the Keplerian energy. 
This relation is displayed in Fig. 5. The value of eccentricity for which the Schwarzschild 
energy parameter is equal to the corresponding Keplerian energy is approximately the same 
as that for which the apocenter distances are equal, eo ~ e+ ~ 1/2. (See Sec. II A.) 
However, this approximate energy parametrization is not expected to be accurate for values 
of e approaching 1/2. A more accurate treatment in Sees. Ill A and IIIB results in e+ < eo, 
as expected. 



III. SELF-CONSISTENT KEPLERIAN LIMIT 



The relativistic equation of orbit (25) is useful for describing characteristics of relativistic 
orbits in a Keplerian limit. However, the solution is inaccurate in describing long-term 
orbital behavior. It is easy to verify that (17) solves (10) only to first-order in e^j. The 
Schwarzschild equation of orbit (17) may instead be expressed as 



— = ^(1 +e cos kip), (44) 

r Tr 



where Tc, e, and k are defined in Eqs. (18)-(21). Substituting (44) into the equation of 
motion (10) and discarding the term nonlinear in eccentricity results in the two conditions: 

^Qej^,{fjrc)-'] (45) 
(fc/rj-i ^ l + 3eA,(fc/rJ-l (46) 

The first condition (45) is satisfied only to first order in ej^j, and the second condition (46) 
is satisfied only to second order in ej^.^. Fortunately, consistency is achieved by taking fc and 
k to be independent functions of e^j and making the following replacements in Eq. (44): 

(47) 

rc/rc ^ rc/rc - \Xl. (48) 
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The constants Kq and Aq are then found by substituting Eq. (44) into the equation of motion 
(10), resulting in 

nl = l- 6e,, - (1 - 126,,)^ = ISe^, + 0(e|,), (49) 

= - (1 - 12%)^ = 544, + 0{et,). (50) 

1 - Se^j 

A self-consistent equation of orbit may now be expressed as 

— = l + ecosfi;a3, (51) 
r 

with the following definitions: 



rVre = | + iVl-12e^,; (52) 



e/e ^ (53) 
k^{l- 126,,)i (54) 



This solution is consistent to all orders in e^, and the more approximate orbital equation 
(25) of Sec. II is recovered by expanding (52)-(54) and keeping terms first order in ej^: 

fe^re(l-3e,,-9e2,-54e|,); (55) 
e^e(l + 3e,,, + 9e|, + 27e^,); (56) 

K ~ 1 _ - _ 189 3 /cryN 

Long-term behavior is predicted very accurately using the exact expressions (52)-(54). The 
radius of relativistic circular orbit (52) is identical to the result obtained by minimizing the 
effective potential (31). 

A systematic verification may be carried out by direct substitution of Eq. (51) into 
Eq. (10) using the exact definitions (52)-(54) for fc, e, and k, and discarding the term 
nonlinear in eccentricity. However, the justification for discarding the term nonlinear in ec- 
centricity is the correspondence principle. Arguments concerning which terms to keep based 
only on a direct comparison of the relative magnitudes of higher-order and lower-order terms 
lead to contradictions. Rather, the domain of validity is expressed by subjecting the solution 
(51) to condition (12). Evaluating the equation of orbit (51) at pericenter results in 

^ = (58) 
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Therefore, according to (12), the domain of vahdity is given by 

e/e + 2(l -fe/rc) < 1, (59) 

Where fc/rc and e are given by (52) and (53), respectively. This condition is consistent with 
Eq. (27) to first order in e^^. Of course, the further condition e^j < 1/12 is also necessary. 

It is worth mention that this self-consistency program is independent of the eccentricity. 
A replacement of the form e ^ e + ?7o, such as in (47) and (48), may be imposed, but rjQ 
cannot be determined by insisting that the solution (44) satisfy the equation of motion (10). 
No adjustments are made to further increase the accuracy of the relativistic eccentricity. 
The accuracy of the relativistic eccentricity is verified only by agreements of (51) with exact 
numerical solutions of the equation of orbit (10). 

The self-consistent equation of orbit (51) is in agreement with exact numerical solutions of 
(10) for a large parameter space, including high-eccentricity (1 > e ^ 0) and very relativistic 
environments (1/12 > e^j ^ 0). These periodic solutions are compared using relative errors 
in radial coordinate and angular frequency separately. The relative error in radial coordinate 
is defined as 

5r/r ^ 'JllZltl^ (60) 
r/r^ 

where r/rc is the exact numerical solution of (10), and f/rc is the self-consistent solution 
(51). For this purpose solutions are compared over a small number or periods. The relative 
error in angular frequency is defined as 

5k/k = 1 - V9n/(^n, (61) 

where Kf/^in = 27rn is the phase of the numerical solution over n cycles, and Kc^n = 27rn is the 
phase of the self-consistent solution over the same number of cycles. This error is determined 
using n = 1600 cycles, and in each case ip^ — (p^ represents the total phase difference. 

Errors for several cases, including relatively large values of e^,^ and e, are plotted in 
Figs. 6-8. For reference: ej^ ~ 10"*^ for the earth-sun system (rc ~ 10®r^j/2); e^^j ~ 10~^ 
for an outlying (rc ~ 1 Mpc) giant spiral galaxy orbiting a cluster containing one thousand 
giant spiral galaxies (rc ~ lO^rjj/2); e^j = 10~^ for a star orbiting a super-massive blackhole 
at a distance corresponding to rc = 500r^,j. It is worth noting that in every case the largest 
error is near apocenter [cp = n), and that this error is in favor of the stated characteristic of 
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reduced apocenter. The exact numerical value for the radial coordinate is smaller than that 
determined by the equation of orbit derived in this Keplerian limit. 

This self-consistent solution (51) is identical in form to the more approximate solution (25) 
of Sec. II. However, the coefficients (52)-(54) result in an equation of orbit that accurately 
predicts long-term behavior. The more approximate solution of Sec. II is a limiting case of 
this self-consistent solution, lending value to the simpler approach and first-order predictions. 



A. Characteristics of Schwarzschild Orbits 



The characteristics of Schwarzschild orbits are expressed more accurately using the self- 
consistent equation of orbit (51)-(54). The self-consistent equation of orbit predicts the rate 
of precession to be 



(27r)~^A(^ = - 1, (62) 

45 .2 
2 



3eM + f4 + (^(4), (63) 



The second-order term Aip^'^^ ~ y^m is in close agreement with that calculated by Ashby 
[21] A^f ~ (f + je^)4i, so that 1.58 < A^(2)/^^(2) ^ ^ q^_ 

The self-consistent equation of orbit predicts the radius of stable circular orbit to be 
reduced, 

SfJr, = -3ej,,-9el,-Oie!,). (64) 

This expression is in agreement with the Schwarzschild geometry; The radius of circular 
orbit as predicted by this self-consistent solution (52) is identical to that calculated from the 
Schwarzschild geometry by minimizing the effective potential (31). The relativistic apsides 
are 

r± r± = — — r, (65) 

r, 1 T e/(l - 3e^j) 

The relativistic apocenter is smaller than the corresponding Keplerian apocenter for values 
of eccentricity less than 

+ (66) 
^1(1 + 14)- (67) 
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(There is no correction first order in e^^ here.) 

The self-consistent equation of orbit predicts the eccentricity of a noncircular orbit to be 
increased, 

5e/e = 3ej^.je/e (68) 
= 3e^, + 9e|, + 0(e|,). (69) 



B. Schwarzschild Energy Parameters 

A more accurate relation between the Schwarzschild energy parameter for a circular orbit 
and that derived from Newtonian mechanics is constructed using the virial theorem together 
with the radius of circular orbit as derived from the self-consistent equation of orbit. (See 
Sec. JIB.) The Schwarzschild energy parameter is given in terms of the Keplerian energy 
by Eqs. (38) and (52), 

= 1 + 2e,, + 9el, + 54e|, + 0{etj). (71) 

This result is also obtained by substituting the radius of circular orbit (31) into the effective 
potential (7), and is in agreement with the approximate relation (39) to first order in e^. 

More generally, a Schwarzschild energy parameter for noncircular, approximately Keple- 
rian orbits is defined by 

e = ;^ ~ — , (72) 

where R± = R±{E) are the relativistic apsides as determined from the intersection of the 
effective potential (7) with a line of constant energy E. That is, 

^''-'"^ -t^-2t- '"W '''' 

The solution is then inverted, yielding E{e). The cubic (73) is solved using Viete's (1540- 
1603) trick [45], 

^ = -^[l + 2p cos + I arccos 7)] , (74) 
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where 0± = vr =1= vr/S, and the following functions are defined: 



f5 = v/1 - 12eM; (75) 
^ ^ P^-^ll _ i8e^^(l - 3ej,,E/E,)]. (76) 



(Notice that 2r^£ ^ = — Therefore, the relativistic eccentricity is given by 

(77) 



\/3/3 sin (| arccos 7) 



1 — /5 cos (| arccos 7) 

This is verified by setting e = and calculating the energy of a circular Schwarzschild orbit, 
resulting in Eq. (70). Inverting (77) results in an expression for the Schwarzschild energy 
parameter, 

E/E^ = 77^ 1 1 - 7^ [1 - cos (3 arccos A)] j , (78) 

where 



3% I ^^^M 



\ = P~\-\l + Vl + ^0, (79) 
C = (80) 

e^^- (81) 

Dividing (78) by (70) and expanding in powers of e^j results in an expression comparable to 
the earlier energy parameterization ansatz (42), 

E/E, = (1 - + 2B'e,, + e\U + 5e')elj ^^^^ 

In the limit e^^ the Keplerian result (41) is recovered. The Schwarzschild energy param- 
eter may instead be expressed in terms of the energy and eccentricity of the corresponding 
Keplerian orbit using (78) together with (53) and (41), 

E/E, = E/E, + 2ejv,(l - 3e^ - e^) ^^^^ 
+ e|,(9 - 18e2 - 37e^ - Se^) + 0(e|,). 

Neglecting terms of orders e^^e^, e|^, and smaller in Eq. (83) results in 

E^E-2ej,,{l-?,e^)\E,\. (84) 
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This result is similar to the energy relation (43) derived from the simple ansatz (42), lending 
value to the much simpler approach of Sec. JIB. Corrections to the orginal ansatz (42) 
are given by Eq. (82). The orginal ansatz (42) and the resulting relation between the 
Schwarzschild energy parameter and corresponding Keplerian energy (43) approximate well 
this more detailed parameterization, (82) and (83), for near-circular orbits requiring small 
relativistic corrections. 

Referring to (83), the energy of a Schwarzschild orbit is smaller than that for the corre- 
sponding Keplerian orbit, until the energies become equal for 

1 - 3e^ - ^ 0, (85) 

so that Co ~ 0.55. The relativistic energy parameter becomes greater than the corresponding 
Keplerian energy for e > Cq. This value of eccentricity is larger than that for which the 
corresponding apocenter distances are equal, as determined from Eq. (67); for e^^ = 1/12, 
e+ ~ 0.52, so that e+ < Cq for all allowed values of ef^j. (Including both the e^^j- and e|^-terms 
in Eq. (83) results in cq ~ 0.550 + 1.99 x 10"^e^j, and the relation e+ < eo holds.) This 
is expected since it is required that f+ > r+ for any energy E = E > Ec, Referring to an 
energy diagram, E must lie below E in order for f+ = r+. 



IV. REISSNER-NORDSTROM ORBITS IN KEPLERIAN LIMIT 

The path of a small, electrically neutral test mass near a spherically-symmetric mass M 
with charge Q is described by the Reissner-Nordstrom geometry [1, 4, 5, 6, 22, 32, 46, 47, 
48, 49, 50]. The metric is given by Eq. (1) with 

where Tq = GQ^/ {AneQC^), and Eq is the electrical permittivity of the vacuum. A far exterior 
solution is described in the Schwarzschild limit, so that r ^ r^j > Tq. Therefore, the 
singularities e^^'-'o^ = are irrelevant in the present context. The equations of motion may 
be expressed as 

i = rV, (87) 
k = {l-r^,/r + r^/r^)i, (88) 
l{e-iy = lr^ + V,s, (89) 
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where an effective potential is defined as 

1 2 2 4- 

2/)-2ta 'c.-L'c C| 'ci 'c /nnX 

Feff = -- + -^-eM^ + eQ^ + ejv,eQ^. (90) 

The parameter Tc is defined in Sec. II after Eq. (7), and the charge-related relativistic 
correction parameter is defined as 



e, 



r2 r2 



(91) 



Q r r 

'm'c 'M 

The condition < r^^ is equivalent to eg < 2e^/, so that for the extreme case in which 
Tq = r^j, Cq = 2ejy,j. However, as discussed below, a consistent treatment of the Reissner- 
Nordstrom (RN) geometry in this Keplerian limit requires the condition eg/ej^j < |. An 
equation for the trajectory of a test particle is obtained by eliminating time from Eq. (89) 
and differentiating once more with respect to if, 

+ ^ = 1 + 3e,, ( - 2e.^ - Ae^.e^ ( '-^t (92) 



Schwarzschild orbits (10) are recovered by setting eg = 0, and the conic-sections of Newto- 
nian mechanics (11) are recovered by setting ej^,^ = eg = 0. 

Reissner- Nordstrom orbits in this Keplerian limit are described by neglecting the last 
terms on the right-hand-sides of Eqs. (90) and (92), resulting in 

1 2 S 2 

iyt I iY>^ itn'-' 

r?r^K,«-^ + i;|-%;i + .g^, (93) 
^''■'+!f«l + 3.„fr^f-2.„!f. (94) 



d<y9^ r r V r / r 

The neglected term is only important when describing very close encounters (r <^ Tc). 
Following the procedure used to solve Eq. (10), described by (12)-(21), a first-order solution 
to (94) may be expressed as 

— ~ 1 -|- e cos (95) 
r 



where 



f,^r,(l-3ejv, + 2eQ), (96) 
e^e(l + 3e,,), (97) 
K^l-3e,, + eQ. (98) 
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A systematic verification may be carried out by substituting (95) into (94), keeping terms 
of orders e, e^j, e^, ee^,j, and ee^ only. However, tlie justification for discarding tlie term 
nonlinear in eccentricity is tlie correspondence principle. Arguments concerning which terms 
to discard based only on direct comparisons of relative magnitudes of higher-order and lower- 
order terms lead to contradictions. Rather, the domain of validity is expressed by subjecting 
the solution (95) to condition (12) at pericenter, resulting in 

e(l + 3e,,) + 2(3e,,-2eQ)«l. (99) 



A. Characteristics of Reissner-Nordstrom Orbits 



Characteristics of Reissner-Nordstrom orbits are described by comparison with Keplerian 
orbits, as in Sec. II A. The approximate equation of orbit (95)-(98) results in the following 
relativistic corrections: 

(27r)-iA<^^3e,,-eQ; (100) 
5fjr, ^ -3ejv, + 26^; (101) 
rVr±^l-3e,,i^ + 2eQ; (102) 

i -)- 6 

5e/e ^ Se^. (103) 

Relativistic corrections due to charge are counter to those due to mass, except that there is 
no charge-related correction to eccentricity. (The exact equation of motion (92) may also 
be linearized, resulting in a charge-related decrease in eccentricity of order 6^/6^.) The rate 
of precession as predicted by the Schwarzschild geometry (29) is decreased when charge is 
present (100). This additional contribution to precession is identical to that calculated by 
Chaliasos [46]. (See also Teli & Palaskar [47].) The radius of circular orbit as predicted by 
the Schwarzschild geometry (30) is increased when charge is present (101). For example, 
in the limit Cg/e^j ^ |, "Tc — ^ Tc- The radius of circular orbit as determined from (101) is 
consistent to first order in e^^j and eg with the stable circular orbit calculated by minimizing 
the effective potential (93), 

^c-Kl + 2eg)r. ^^^^^ 
+ i(l + 26g)r,[l-12e,,(l + 2eg)-2]i 
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The RN effective potential (93) is compared to that derived from Newtonian mechanics (9) 
in Fig. 9. (The Schwarzschild effective potential is included for reference.) The relativistic 
pericenter is reduced for ^ < ^qI^u < §• The remaining parameter space is considered in 
two cases: (Case I) For < Cg/e^/ < \, the relativistic pericenter is enlarged only for very 
large values of eccentricity, | < e_ < 1, for which the approximate equation of orbit (95) 
and resulting orbital characteristics (100)-(103) are not expected to be accurate; (Case II) 
For large values | < Cg/eA/ < 2 it is necessary to include the higher-order ej^/jeg-term in 
the equation of orbit (92) in order to make a consistent argument concerning the value 
of eccentricity beyond which the relativistic pericenter is enlarged when compared to the 
corresponding Keplerian pericenter. This is avoided by defining a Keplerian limit including 
the condition Cg/e^^ ^ 1, or equivalently (jg/r^^)^ ^ 1/2. (It is only necessary to impose 
that Cg/eA/ < |. However, it is convenient to restrict the problem to a smaller parameter 
space because the RN (93) and Newtonian (9) effective potentials intersect at r = ri given 
by Tc/vi = eg/e^^. The problem is simplified if rc/n ^ 1, in which case ri ^ r for all r 
consistent with the domain of validity as given by (12).) Having addressed both cases, it 
may be stated that in this Keplerian limit the relativistic pericenter is always reduced. The 
relativistic apocenter is smaller than the corresponding Keplerian apocenter for eccentricities 
smaller than e+ ~ | — C(eg/eA^). 



B. Reissner-Nordstrom Energy Parameters 

For a particular Keplerian orbit, identified by total energy, it is necessary to identify the 
corresponding Reissner-Nordstrom energy for an orbit described by the same angular mo- 
mentum. This provides a relativistic correction to Newtonian energies for Keplerian orbits. 
This is also useful when comparing orbital properties using energy diagrams. The gener- 
alized virial theorem [39, 40, 41, 42] provides this relation for circular Reissner-Nordstrom 
orbits. Referring to (89) and (90), a RN potential energy parameter for a circular orbit is 
defined as 

GM GMrl GMr, , 
= - %^ + eg^^. (105) 



According to the virial theorem, 

GM GMr^ ^ GMr, 



Tc = ^ + - eg-z^, (106) 
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GMr,( 




2rc fc V 


- f — 

HI ~2 




-2eQ), 



where Tc is a RN kinetic energy parameter. Therefore, using Ec = Tc + Vc, 

GM GMrl , , 

^c = -^ + eM^^ (107) 

(108) 

(109) 

where fc ~ rc(l — 3e^^ + 2eQ) is used in the last step. An approximate RN energy parameter 
is defined by ansatz as in (42), resuhing in a relation between noncircular RN and Keplerian 
energies, 

E^E- [2e,,(l - 4e2) - 26^(1 - e^)]!^^. (HO) 

(This ansatz is investigated in Sec. IV C 2, wherein Eq. (110) is shown to approximate well a 
more detailed parameterization and is a good approximation for near-circular orbits.) The 
RN energy parameter is smaller than the corresponding Keplerian energy for eccentricities 
smaller than cq ~ | + O^eq/ej^j). The value of eccentricity for which the RN energy param- 
eter is equal to the corresponding Keplerian energy is approximately the same as that for 
which the apocenter distances are equal, cq ~ e+ ~ 1/2. (See Sec. IV A.) However, this ap- 
proximate energy parametrization is not expected to be accurate for values of e approaching 
1/2. A more accurate treatment in Sees. IV CI and IV C 2 results in e+ < cq, as expected. 



C. Self- Consistent Keplerian Limit 

A more accurate self-consistent equation of orbit is derived for the Reissner-Nordstrom 
geometry by following the procedure outlined in Sec. Ill, Eqs. (44)-(50). The resulting 
equation of motion may be expressed as 



— = 1 + e cos K(p, 
r 



with the following definitions: 



1 1 

2 + 2 



e/e 



12e, 



'M 



'l + 2er 



1 - 3% ' 



K = {l + 2er 



12e 



(1 + 26. 



;iii) 

;ii2) 
;ii3) 

;ii4) 
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This solution is consistent to all orders in and eg, and the more approximate orbital 
equation (95) is recovered from (111) by expanding (112)-(114) and keeping terms first 
order in both e^^ and eg: 

^ rc(l - 3e^j + 2eQ + Ge^jeg - 9e|^); (115) 
e^e(l + 3e,, + 9e2,); (116) 
k^l- Se^ + eg + Q^Aieg - y4i - (117) 

(There is no correction of order eg in the approximate radius of circular orbit fc.) A sys- 
tematic verification may be carried out by direct substitution of Eq. (Ill) into Eq. (94) 
using the exact definitions (112)-(114) for fc, e, and k, and discarding the term nonlinear 
in eccentricity. However, the justification for discarding the term nonlinear in eccentricity is 
the correspondence principle. Arguments concerning which terms to keep based only on a 
direct comparison of the relative magnitudes of higher-order and lower-order terms lead to 
contradictions. Rather, the domain of validity is expressed by subjecting the solution (111) 
to condition (12). Evaluating the equation of orbit (111) at pericenter results in 

Vc I + e 



r_ rc/r^ ' 

Therefore, according to (12), the domain of validity is given by 



(118) 



e/e + 2(l -fe/rc) < 1, (119) 

which is consistent with Eq. (99) to first order in e^j and eg. This condition reduces to that 
for the self-consistent Schwarzschild solution (59) when eg = 0. 

1. Characteristics of Reissner- Nordstrom Orbits 

The characteristics of Reissner- Nordstrom orbits are expressed more accurately using the 
self-consistent equation of orbit (111)-(114). The self-consistent equation of orbit (111) 
predicts the rate of precession (28) to be 

{2-KY^^^ = ?,eM- eg-lhe„[eg ^ 

The effect of charge is to reduce the rate of precession predicted by the Schwarzschild 
geometry. 
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The self-consistent equation of orbit predicts the radius of stable circular orbit to be 
reduced by 

Srjr, = -3e,, + 26^ + Ge^.e^ - Qe^, + 0{c-'). (121) 

(There is no correction of order eq here.) This expression is in agreement with the Reissner- 
Nordstrom geometry; The radius of circular orbit, as predicted by this self-consistent solution 
(112), is identical to that calculated from the RN geometry (104) by minimizing the effective 
potential (93). The relativistic apsides are 

fjr^ = !jL L£f ^. (122) 

The relativistic apocenter is smaller than the corresponding Keplerian apocenter for values 
of eccentricity less than 

e+ = i[l - KeQ/e^) - 1(^0/%)' ^^^^^ 
+ FM-H + 0{eyel,) + 0{c~% 

(There are no corrections first order in e^,^ and eg here.) 

There is no charge- related correction to eccentricity; 5e/e is identical to that derived from 
the Schwarzschild geometry, Eqs. (68) and (69). 



2. Reissner- Nordstrom Energy Parameters 

A more accurate relation between the Reissner-Nordstrom energy parameter for a circular 
orbit and that derived from Newtonian mechanics is constructed using the virial theorem 
together with the radius of circular orbit as derived from the self-consistent equation of 
orbit. (See Sec. IV B.) The Reissner-Nordstrom energy parameter is given in terms of the 
Keplerian energy by Eqs. (108) and (112), 

(l + 26Q)(rc/rc)-4eM 



Ec/Ec 



(124) 



1 + 2ej,^ - 2eQ - 12eMeQ 



+ 9e^, + 46^ + 0(0. 

This result is also obtained by substituting the radius of circular orbit (112) into the effective 
potential (93). This expression is in agreement with the approximate relation (109) to first 



order in e^j and eg. 
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More generally, a Reissner-Nordstrom energy parameter for noncircular, approximately 
Keplerian orbits is defined by 

— i?„ f j„ — f_ 



R+ + R_ f , + r 



(125) 



where R± = R±{E) are the relativistic apsides as determined from the intersection of the 
effective potential (93) with a line of constant energy E. That is, 

^ ly-t 1 ,*-j2 ,j^3 lyi^ 

The solution is then inverted to give E{e). The cubic (126) is solved using Viete's (1540- 
1603) trick [45], 

r 1 + 2e 

^ = — ^ — ^[l + 2/?cos(0± + iarccos7)l, (127) 
R± 6% 

where = tt =l= 7r/3, and the following functions are defined: 



[l-126,,(l + 2eQ)-2]^ (128) 



18% Se^^ E 



(1 + 26^)2 V" l + 2eaE, 
(Notice that 2r'^i~'^ = —l/Ec-) Therefore, the eccentricity is given by 

Vsp sin ( I arccos 7) 



(129) 



(130) 



1 — (3 cos (3 arccos 7) 

This is verified by setting e = and calculating the energy of a circular RN orbit, resulting 
in Eq. (124). Inverting (130) results in an expression for the RN energy parameter 



3% I ISe^f ^131) 



X [1-/9^ cos (3 arccos A)] 

where 



A = r'r'(i + v^T+ce), (132) 

C = (133) 



e^^. (134, 
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Dividing (131) by (124) and expanding in powers of e^y^j and eg results in an expression 
comparable to the earlier energy parameterization ansatz (42), 

E/E, = (l-e~2)[l + 2e~\, 

(135) 

+ 6^(14 + he^)el, - ^e\,,eQ\ + 0{c-'). 

(There is no term of order eg here.) The Schwarzschild result (82) is recovered by setting 
eg = 0, and the Keplerian result (41) is recovered by setting e^^ = eg = 0. The RN 
energy parameter may instead be expressed in terms of the energy and eccentricity of the 
corresponding Keplerian orbit using (131) together with (113) and (41), 

E/E, = E/E, + 2e^,(l - 36^ - e^) - 2eg(l - e^) 



„2 

(136) 

6^2,(9- 18e2-37e^-5e6) 



Neglecting terms of orders e^^e^, e^^eg, e|^, eg, and smaller in Eq. (136) results in 

E^E- [2ejy,,(l - 36=^) - 2eg(l - e^)]\E,\. (137) 

This result is similar to the energy relation (110) derived from the simple ansatz (42), 
lending value to the much simpler approach of Sec. IV B. Corrections to the orginal ansatz 
(42) are given by Eq. (135). The orginal ansatz (42) and the resulting relation between 
the RN energy parameter and corresponding Keplerian energy (110) approximate well this 
more detailed parameterization, (135) and (136), for near-circular orbits requiring small 
relativistic corrections. 

Reffering to (136), the energy of a Reissner- Nordstrom orbit is smaller than that for the 
corresponding Keplerian orbit, until the energies become equal for 

26,,(1 - 3eg - et) - 26g(l - el) ^ 0, (138) 

so that cq ^ 0.55 — 0.44(eg/ej,f). The relativistic energy parameter becomes greater than 
the corresponding Keplerian energy for e > cq. For ^q/^m ^ 1 ^^is value of eccentricity 
(eo) is larger than that for which the corresponding apocenter distances are equal (e+), as 
determined from Eq. (123). For example, choosing eg/% = 1/10 results in e+ ~ 0.48 and 
Co ~ 0.51. This is expected since it is required that f+ > for any energy E = E > Ec] 
Referring to an energy diagram, E must lie below E in order for f_|_ = r_|_. 
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V. SCHWARZSCHILD-DE SITTER ORBITS IN KEPLERIAN LIMIT 



The path of a small test mass near a spherically-symmetric central mass M including 
the effect of the cosmological constant A is described by the Schwarzschild-de Sitter (SdS) 
geometry [3, 4, 5, 20, 32, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65]. The metric 
is given by Eq. (1) with 

where = 3/A, and A > (repulsive). A solution is described in a Keplerian limit, so that 
Tf^^ r ^ r]^,j. The radius of the cosmological horizon is large when compared to the scale 
of the orbits. Therefore, the singularities e^^^^o) = Q are irrelevant in the present context. 
The equations of motion may be expressed as 

i = r^Lf, (140) 
k = {l-r,Jr-ryrl)i, (141) 
i(fc2-l)c2 = if2 + V;fr, (142) 

where an effective potential is defined as 

1 2 S 2 1 2 

ly I itrt^ ly'-^ ly'-' I ly"-' 

r^r^K,^ -f + 5# - -Mf, - 'a;5 - ("3) 

The A-related relativistic correction parameter is defined as 

- 772 = ¥Mf,- (144) 
'm'k 'a 

It is assumed that <^ e^^, or equivalently A <^ r^j/r^. This is reasonable considering the 
very small value for the cosmological constant [66, 67, 68]: A ^ Hq/c^ «i 10~^^cm~^. An 
equation for the trajectory of a test particle is obtained by eliminating time from Eq. (142) 
and differentiating once more with respect to (p, 

+:i^l+3,X-^-2eX-r. (146) 



d<^^ r r \ r / V r 

Schwarzschild orbits (10) are recovered by setting = 0, and the conic-sections of Newtonian 
mechanics (11) are recovered by setting e^;.^ = = 0. 

Although a self-consistent solution for SdS orbits in this Keplerian limit is intractable, 
a first-order solution provides three orbital characteristics with much less effort. Following 

26 



the procedure used to solve Eq. (10), described by (12)-(21), a first-order solution to (145) 
may be expressed as 



1 + ecos Koj, (146) 
r 



where 



fc^r,(l-3eA, + 2eJ, (147) 
e^e(l + 3e^, + 8ej, (148) 
kp^l- Se^j - Se^. (149) 

A systematic verification may be carried out by substituting (146) into (145), keeping terms 
of orders e, e^^^, e^, ee^^^, and ee^ only. However, the justification for discarding the term 
nonlinear in eccentricity is the correspondence principle. Arguments concerning which terms 
to discard based only on direct comparisons of relative magnitudes of higher-order and lower- 
order terms lead to contradictions. Rather, the domain of validity is expressed by subjecting 
the solution (146) to condition (12) at pericenter, resulting in 

e(l + Se^ + 8eJ + 2(36^ - 26^) < 1. (150) 
A. Characteristics of Schwarzschild-de Sitter Orbits 

Characteristics of Schwarzschild-de Sitter orbits are described by comparison with Kep- 
lerian orbits, as in Sec. II A. The approximate equation of orbit (146)-(149) results in the 
following relativistic corrections: 

(27r)-iAv^^3ejv, + 3eA; (151) 

5fjr, ^ -3e^, + 26^; (152) 

1^2e l±3e , , 

r±/r± ^ 1 - 36,,^^ + 2eAy^; (153) 

5e/e ^ 3ej,j + 8e^. (154) 

Relativistic corrections due to A include both increased rate of precession and increased 
eccentricity. (Recall from Sec. IV A that there is no correction to eccentricity due to charge.) 
The rate of precession as predicted by the Schwarzschild geometry (29) is increased if A is 
present (151). The additional contribution, (27r)^^Ay9^ ^ 3e^, is in agreement with the 
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standard perturbative result [3, 63]. The eccentricity as predicted by the Schwarzschild 
geometry (34) is increased if A is present (154). Another effect of A is to increase the overall 
size of the orbit, when compared to Schwarzschild orbits. The radius of circular orbit as 
predicted by the Schwarzschild geometry (30) is increased when A is present (152). The 
present formalism provides a determination of the radius of relativistic circular orbit to first 
order in and (147); the standard approach of minimizing the effective potential results 
in a quintic equation for fc. 

The relativistic pericenter (153) is always reduced. The condition for which f_ = r_ 
is e^/e^/ > 3/2, which is inconsistent with the stated assumption of this Keplerian limit: 
^ When compared to a Schwarzschild orbit the effect of A is to increase the 
pericenter for values e < 1/3 and to decrease the pericenter for values e > 1/3. The 
relativistic apocenter is smaller than the corresponding Keplerian apocenter for eccentricities 
smaller than e+ ~ 1 — C(e^/e^j). 



B. Schwarzschild-de Sitter Energy Parameters 

For a particular Keplerian orbit, identified by total energy, it is necessary to identify 
the corresponding Schwarzschild-de Sitter energy for an orbit described by the same angular 
momentum. This provides a relativistic correction to Newtonian energies for Keplerian 
orbits. This is also useful when comparing orbital properties using energy diagrams. The 
generalized virial theorem [39, 40, 41, 42] provides this relation for circular Schwarzschild- 
de Sitter orbits. Referring to (142) and (143), a SdS potential energy parameter for a circular 
orbit is defined as 

GM GMrl GMfl 
Vc = — - Cm^^ - Ca (155) 



' c ' c 



According to the virial theorem. 



GM GMr^ GMf"^^ , , 

Te = + 3e,,— 3^ - eA^^, (156) 
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where Tc is a SdS kinetic energy parameter. Therefore, using Ec = Tc + Vc, 

~ GM GMrl ^ GMfl 

^c = -^ + %.^-2e,^ (157) 

a £,(1 + 2e„ + 2ej), (169) 

where fc ~ rc(l — 3eJ^/ + 2€j^) is used in the last step. An approximate SdS energy parameter 
is defined by ansatz as in (42), resuhing in a relation between noncircular SdS and Keplerian 
energies, 

E^E- [2eA,(l - 4e2) + 26^(1 - 9e^)]\E,\. (160) 

The SdS energy parameter is smaller than the corresponding Keplerian energy for eccentric- 
ities smaller than eo ~ | — 0{ejy/ej^.j). The value of eccentricity for which the SdS energy 
parameter is equal to the corresponding Keplerian energy is approximately the same as 
that for which the apocenter distances are equal, eo ~ e+ ~ 1/2. (See Sec. VA.) How- 
ever, this approximate energy parametrization is not expected to be accurate for values of 
e approaching 1/2. 



VI. SUMMARY 



The relativistic central-mass problem is investigated in a Keplerian limit. Beginning 
with the Schwarzschild metric, a relativistic equation of orbit is derived that is similar in 
form to that describing Keplerian orbits. This equation of orbit includes three relativistic 
corrections to Keplerian orbits: precession; reduced radial coordinate; and increased eccen- 
tricity. The prediction for the relativistic contribution to precession is in agreement with 
existing calculations and observations. The predicted reduction in size of a circular orbit is 
also in agreement with existing calculations. These agreements provide confidence in a new 
quantitative prediction of increased eccentricity, which may be subjected to observational 
tests. 

The methods and approximations describing this Keplerian limit to the Schwarzschild 
geometry are also applied to the Reissner- Nordstrom and Schwarzschild-de Sitter metrics. 
The resulting equations of orbit are identical in form to that derived for the Schwarzschild 
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metric and include relativistic corrections due to charge (RN) and cosmological constant 
(SdS). In every case the relativistic equation of orbit has a form that is easily compared to 
that describing Keplerian orbits (11) of Newtonian mechanics, 

T 

— = l + ecosKv^. (161) 
r 

The coefficients fc, e, and k provide relativistic corrections for each geometry. These cor- 
rections are given to first order in Table I. The first-order shift in apside (precession) per 

TABLE I: Summary of first-order relativistic corrections to Keplerian (K) orbits for the 
Schwarzschild (S), Reissner-Nordstrom (RN), and Schwarzschild-de Sitter (SdS) geometries. The 
coefficients of the relativistic equation of orbit, Eq. (161), are given for each of the three geometries. 
In the bottom two rows the subscript X represents mass M, charge Q, or cosmological constant A. 
The parameter r^. = jGM is the radius of a circular Keplerian orbit with angular momentum I. 
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1 - 3ej^^ 


1 - 3ej.; + eg 


1 - 3e,, - 3e^ 


fc/rc 
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1 - 3ej,^ 


1 - 36^^ + 2eQ 


1 - 3e;v, + 2e^ 


e/e 


1 


l + 3e,, 


1 + 36,, 


1 + 36,, + 86a 
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rq 


_ '''c 
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2GM 


''^ ~ Uvreoc^^ 





revolution A(/? = 27r(K~^ — 1) is consistent with known results for the Schwarzschild, Reissner- 
Nordstrom, and Schwarzschild-de Sitter geometries. (See Table II, left column.) The RN 
geometry predicts an additional contribution to precession in the opposite direction to that 
due to matter, while the SdS geometry predicts an additional contribution to precession in 
the same direction as that due to matter. For each of the three geometries, the radius of 
circular orbit is consistent to first order with that determined by minimizing the relativistic 
effective potential. The Schwarzschild geometry predicts a reduced radius of circular orbit. 
Both the RN and SdS geometries predict an radius of circular orbit that is larger than that 
predicted by the Schwarzschild geometry, but still smaller than that for a corresponding 
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Keplerian orbit. Schwarzschild orbits are predicted to be more eccentric than corresponding 
Keplerian orbits, and a cosmological constant (SdS) serves to further increase the eccentric- 
ity. The presence of electric charge (RN) does not result in any additional contribution to 
eccentricity. Relativistic corrections to eccentricity may serve as additional tests of general 
relativity. 

TABLE II: Summary of additional first-order relativistic corrections to Keplerian (K) orbits for 
the Schwarzschild (S), Reissner-Nordstrom (RN), and Schwarzschild-de Sitter (SdS) geometries. 
First-order angular shifts in apsides (precession) per revolution are listed in the left column. First- 
order relativistic corrections to apsides are listed in the right column. The relativistic correction 
parameters are listed in Table I, bottom two rows. 
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This model and the resulting first-order corrections are valid for near-circular orbits 
(e ^ 1/2) that require only small relativistic corrections (e^^ <^ 1/12; eg ^ ^ ^m)- 

In addition to the properties listed in Table I, the overall size of a Schwarzschild orbit is 
predicted to be smaller than a corresponding Keplerian orbit. This is determined not only 
by the radius of circular orbit, but also by a comparision of relativistic apsides to those for 
corresponding Keplerian orbits. Both the apocenter and pericenter distances are found to 
be smaller for Schwarzschild orbits. (See Table II, right column.) Both the RN and SdS 
geometries predict apsides that are larger than the Schwarzschild apsides, but still smaller 
than the corresponding Keplerian apsides. 

Long-term orbital behavior is predicted very accurately using a self-consistent Keplerian 
limit. This is demonstrated by comparing a self-consistent equation of orbit to the exact 
numerical solution for Schwarzschild orbits. The self-consistent equation of orbit is also 
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given by Eq. (161), but with more accurate expressions for the coefficients fc, e, and k. For 
examples, the self-consistent solution predicts: a radius of circular orbit that is identical 
to that calculated by minimizing the relativistic effective potential; and a relative error in 
angular frequency of 10~^° over 1600 cycles. (See Fig. 6.) This solution is accurate in 
describing Schwarzschild orbits with eccentricities as large as e = 1/2 and requiring large 
relativistic corrections. The coefficients in Table I are found to be limiting cases of those 
derived for the more accurate self-consistent equation of orbit, lending value to the simpler 
approach. A self-consistent equation of orbit is also derived for the Reissner-Nordstrom ge- 
ometry, foregoing numerical studies. A self-consistent model for the Schwarzschild-de Sitter 
geometry is intractable and is not pursued. Very accurate relativistic energy parameters for 

TABLE III: Summary of relations between Newtonian energies for Keplerian (K) orbits and 
Schwarzschild (S), Reissner-Nordstrom (RN), and Schwarzschild-de Sitter (SdS) energy parame- 
ters. First-order relativistic corrections to energies for circular orbits are listed in the left column. 
First-order relativistic corrections to energies for near-Keplerian orbits are listed in the right col- 
umn. The relativistic correction parameters are listed in Table I, bottom two rows. 
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circular orbits are derived using the virial theorem. This is useful for comparing the energy 
of a circular relativistic orbit Ec to the energy of a corresponding circular Keplerian orbit 
Ec- This relation is summarized for the three geometries in Table III, left column. (See 
also Fig. 4.) Because the relativistic orbits are taken to be very near-Keplerian, an energy 
parameterization for noncircular bound orbits analogous to that for Newtonian mechanics 
is investigated. The total energy of a relativistic orbit E is defined by simple ansatz, 

E = {1- e^)E„ (162) 
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where Ec is the energy of a relativistic circular orbit, and e is a new relativistic eccentricity 
derived in the context of the relativistic equation of orbit (161). Then, using the relation 
between e and e (Table I, middle row), together with the relation between Ec and Ec 
(Table III, left column), a relation between the total energy for a noncircular relativistic 
orbit and total energy for a Keplerian orbit is derived. This relation is summarized for 
the three geometries in Table III, right column. (See also Fig. 5.) The virial theroem 
and simple ansatz (162) provide first-order relativistic corrections to Newtonian energies for 
bound orbits. Finally, this simple energy parameterization is compared to a more detailed 
parameterization constructed using the intersection of the relativistic effective potential 
with a line of constant energy. The results are similar, lending value to the simple ansatz 
(162) and resulting approximate relations. A more detailed energy parameterization for the 
Schwarzschild-de Sitter geometry results in a quintic equation for E and is not pursued. 

The additional unstable circular orbit and relativistic capture described in the standard 
Newtonian limit to general relativity [1, 3, 5, 11, 22, 32, 35] are absent in the present treat- 
ment. However, the present approach to the relativistic central- mass problem results in an 
equation of orbit that exhibits several characteristics of relativistic orbits at once. In this 
Keplerian limit characteristics of general-relativistic orbits are provided as corrections to 
Keplerian orbits of Newtonian mechanics, providing a qualitative and quantitative under- 
standing of the effects of general relativity on bound systems. It should also be possible 
to adapt these results to relative motion of binary systems [56, 61]. Perhaps more gen- 
eral statements may be made concerning larger systems and more extreme environments. 
Globular clusters are expected to be biased toward high eccentricity in galaxies with larger 
cores. Individual stars orbiting near blackholes are expected to be in anomolously small and 
eccentric orbits. These characteristics could also serve as an indicator of dark matter and 
dark energy. Consider two galaxies, each having approximately the same amount of visible 
matter. A large difference in the amount of coexisting dark matter should be apparent in the 
eccentricities of the orbits of individual stars and star clusters. Individual outlying members 
of galaxy clusters are expected to be biased toward high eccentricity and large precession 
rates due to both the large central mass and the cosmological constant. Although effects of 
the cosmological constant on planetary orbits have been ruled out, they may be observable 
in galaxy clusters with larger radii, for which ej^ ~ Ar^ becomes non-negligible. 

The methods and approximations describing this Keplerian limit may be applied to other 



33 



static spherical spacetimes. The results summarized in Tables I, II and III already describe 
related geometries using simple replacements. The results for the Reissner-Nordstrom geom- 
etry are extended to include magnetic charge with the replacement Q^/^o ~^ Q^/^o + -P^/^oC^) 
where P is the magnetic charge [4, 22]. The results for the Schwarzschild-de Sitter geometry 
become those for the Schwarzschild-anti-de Sitter (A < 0) with the replacement e^^ — 
[23, 52, 54, 69]. It may also be possible to apply this Keplerian limit to more exotic objects 
such as wormholes [70], naked singularities, and Boson and Fermion stars [71]. 
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FIG. 1: Central- mass orbit derived from the Schwarzschild geometry (solid), Eq. (25), compared 
to a corresponding Keplerian orbit (dashed), Eq. (11). Precession is one orbital characteristic due 
to relativity, and is illustrated here for < (p < 20tt. The eccentricity is chosen to be e = 0.2 for 
both the Schwarzschild and Keplerian orbits. This effect is exaggerated by the choice of relativistic 
correction parameter (e^j- = 0.1) for purposes of illustration. However, precession is present for 
smaller (non-zero), reasonably chosen values of e^j as well. 
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FIG. 2: Effective potential (7) derived from the Schwarzschild geometry (solid) compared to that 
derived from Newtonian mechanics (dashed). The vertical dotted lines identify the radius of circular 
orbit as predicted by Schwarzschild Rc, and Newton rc- The Schwarzschild geometry predicts a 
smaller radius of circular orbit as given by Eq. (31). The value = 0.06 has been chosen 
for purposes of illustration. This large value is inconsistent with the condition e^^ <^ 1/12 and, 
therefore, with the approximation of Eq. (30). However, a self-consistent solution derived in Sec. Ill 
correctly predicts the radius of circular Schwarzschild orbit (52) for this large relativistic correction 
parameter. 
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FIG. 3: Comparison of Schwarzschild orbit (25) to a corresponding Keplerian orbit (11). Reduced 
radial coordinate is one orbital characteristic due to relativity. The semi-minor axis is reduced 
more than the semi-major axis, so that the Schwarzschild orbit is also more eccentric than the 
corresponding Keplerian orbit. These effects are exaggerated by the choice of parameters (e = 
0.4, €j^,j = 0.1) for purposes of illustration, and precession has been removed (k 1) in order to 
emphasize the size and shape of the Schwarzschild orbit. However, these effects are present for 
smaller (non-zero), reasonably chosen values of e and e^j as well. 
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FIG. 4: Effective potential (7) derived from the Schwarzschild geometry (solid) compared to that 
derived from Newtonian mechanics (dashed). Orbital energies are superimposed using dotted 
horizontal lines; The energy for a circular orbit as predicted by Newtonian mechanics (top) is 
greater than the Schwarzschild energy parameter (bottom) describing the corresponding relativistic 
circular orbit. The energies are related approximately by Eq. (40). The value ej^j = 10~^ has been 
chosen. 
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FIG. 5: Effective potential (7) derived from the Schwarzschild geometry (solid) compared to that 
derived from Newtonian mechanics (dashed). Orbital energies are superimposed using dotted 
horizontal lines; For the chosen value of eccentricity, the energy for an elliptical orbit as predicted by 
Newtonian mechanics (top) is greater than the Schwarzschild energy parameter (bottom) describing 
the corresponding noncircular Schwarzschild orbit. The energies are related approximately by 
Eq. (43). The values e = 0.2 and e^^ = 10~^ have been chosen. 
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FIG. 6: Comparison of self-consistent Schwarzschild solution (51) to the exact numerical solution 
of (10) for e^^ = 10~^ and moderate values of e. The relative error, Eq. (60), is plotted for the 
initial interval < (/9 < 47r. The phases are nearly equal over this initial small interval, so that this 
figure effectively represents the relative error in the radial coordinate: 5r/r ~ 10^^ for e = 0.1; 
and 5r/r^ 10~^ for e = 0.2. In each of these two cases the relative error in angular frequency, 
Eq. (61), is found to be: Sk/k ~ 10"^°. 
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FIG. 7: Comparison of self-consistent Schwarzschild solution (51) to the exact numerical solution 
of (10) for intermediate values of e^^ and e. The relative error, Eq. (60), is plotted for the initial 
interval < (p < Att. The phases are nearly equal over this initial small interval, so that this figure 
effectively represents the relative error in the radial coordinate: 5r /r ~ 10~^. In each of these two 
cases the relative error in angular frequency, Eq. (61), is found to be: ~ 10~^. 
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FIG. 8: Comparison of self-consistent Schwarzschild solution (51) to the exact numerical solution 
of (10) for more extreme values of ej^/j and e. The relative error, Eq. (60), is plotted for the initial 
interval < (f < 47r. The phases are nearly equal over this initial small interval, so that this figure 
effectively represents the relative error in the radial coordinate: (5r/r ~ 10"^. The relative error 
in angular frequency, Eq. (61), is, for each of these two cases: 6k/k ~ 10^^ for e^j = 10~^ and 
e = 0.8; and 5k/k ~ 10^*^ for ej^,j = 10"^ and e = 0.4. 
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FIG. 9: Effective potential (93) as derived from the Reissner-Nordstrom geometry (solid) compared 
to that derived from Newtonian mechanics (upper, long dashes) and Schwarzschild geometry (lower, 
short dashes). The vertical dotted lines identify the radius of circular orbit as predicted by (from 
left to right) Schwarzschild, Reissner-Nordstrom Rc, and Newton rc- For the chosen values of rela- 
tivistic correction parameters, the Reissner-Nordstrom geometry predicts a radius of circular orbit 
(104) that is smaller than the corresponding Newtonian orbit and larger than the corresponding 
Schwarzschild orbit. The values e^^j = 0.06 and eg = 0.02, for which {^q / rj^j)^ = 1/6, have been cho- 
sen for purposes of illustration. These large values are inconsistent with the condition eg/e^,/ <C 1 
and, therefore, inconsistent with the approximation of Eq. (101). However, a self-consistent solu- 
tion derived in Sec. IV C correctly predicts the radius of circular Reissner-Nordstrom orbit (112) 
for these large relativistic correction parameters. 
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